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Abstract. We derive a posteriori error bounds for a quasilincar parabolic 
problem, which is approximated by the hp- version interior penalty discontinu- 
ous Galerkin method (IPDG). The error is measured in the energy norm. The 
theory is developed for the semidiscrete case for simplicity, allowing to focus 
on the challenges of a posteriori error control of IPDG space-discretizations 
of strictly monotone quasilinear parabolic problems. The a posteriori bounds 
are derived using the elliptic reconstruction framework, utilizing available a 
posteriori error bounds for the corresponding steady-state elliptic problem. 



Discontinuous Galerkin (DG) methods |Bak77| IWhe78| IArn82] . have enjoyed 
substantial development in recent years. For parabolic problems DG methods are 
interesting due to their good local conservation properties as well as due to their 
block-diagonal mass matrices. 

This work is concerned with the derivation of a posteriori error bounds for the 
space-discrete interior penalty discontinuous Galerkin method (IPDG) for quasilin- 
ear parabolic problems with strictly monotone non-linearities of Lipschitz growth. 

A posteriori error bounds for ft,- version DG methods are derived in [KP03, KP07, 
IBHL031 IHSW07] and for DG-in-space parabolic problems in [SW051 IEP051 IYC061 
IGL08 , EV09. The contribution of this work is twofold: 




• the derivation of a posteriori energy-norm error bounds for IPDG methods 
for quasilinear parabolic problems, and 

• the resulting a posteriori bounds are are explicit with respect to the local 
elemental polynomial degree. 

A key tool in our a posteriori error analysis is the elliptic reconstruction technique 
|MN03| l"LM06| IGL08] . Roughly speaking, in the elliptic reconstruction framework 
the error is split into a parabolic and an elliptic part, respectively. In the interest 
of being explicit with respect to the dependence of the a posteriori error bounds 
in the elemental polynomial degree p, we restrict the presentation to quadrilateral 
elements of tensor-product type (cf. Remark 13.71) . 



Let f2 be a bounded open (curvilinear) polygonal domain with Lipschitz bound- 
ary dfl in R d , d = 2, 3. For u C SI, we consider the standard spaces L 2 {uj) (whose 
norm is denoted by || • for brevity), H 1 ^) and Hq(oj), whose norm will be denoted 
by || ■ ||i, along with its dual with norm ||-||_i. For brevity, the standard 

inner product on L 2 (f2) will be denoted by (•, •) and the corresponding norm by || -|| . 



1. Introduction 




2. Model problem and the IPDG method 
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We also define the spaces L 2 {0,T,X), X £ {L 2 {uj), H ±l {ui)} and L°°(0, T, L 2 (51)), 
consisting of all measurable functions v : [0,T] — > X, for which ||u||z/ 2 (o,T , ;jf ) : — 

(Jo ll v Wllx) 1/2 < +°° and |MU~(o,T ; L2(fi)) := ess sup te[0jT] \\v{t)\\. (The differ- 
entials in the integrals with respect to t are suppressed for brevity throughout this 
work.) 

We identify function v £ [0, T] x 51 — > R with w : i — > X and we denote u(t), 
t G [Q,T], for »e[0,r]xft^K. 

For t € (0,T], we consider the problem of finding a function it satisfying 

(1) u t (t,x) - V- (a(t,x, \Vu(t,x)\)Vu(t,x)) = f(t,x) in (0,T] x 51, 

where / G L° C (Q, T; L 2 (51)) and a scalar uniformly continuous function, subject to 
initial condition u(0,x) = uq(x) on {0} x 51, for uq G L 2 (51), and homogeneous 
Dirichlet boundary conditions on [0, T] x 951. 

We assume that the non-linearity a in equation ([TJ is of strongly monotone type 
with Lipschitz growth so that there exist positive constants a and a such that the 
following inequalities hold: 

(2) \a(t,x,\y\)y - a(t,x,\z\)z\ < a\y - z\ 

(3) (a(t,x, \y\)y - a(t,x, \z\)z) ■ (y - z) > a\y-z\ 2 , 

for all vectors y, z G R d , and all (i, x) G [0, T] x 51. 

Let T be a shape-regular subdivision of 51 into disjoint closed quadrilateral el- 
ements k £ T. We assume that n £ T are constructed via C M -diffeomorphisms 
with non-singular Jacobian F K : (—1, 1) — > k, so as to ensure 51 = U Ke T-K. 

For p £ N, Q p (k) is the set of all tensor-product polynomials on (— 1, l) d of 
degree p in each variable and let 

(4) S p := {v £ L 2 (S1) : v\ Fk £ Q p {k), k £ T}, 

be the (discontinuous) finite element space. Let T be the union of all (d — 1)- 
dimensional element faces e associated with the subdivision T (including the bound- 
ary). Let also r int := r\951, so that T = 951 U r int . 

Let k + , k~ be two (generic) elements sharing a face e := k + l~l kT C Tint with 
respective outward normal unit vectors n + and n on e. For q : 51 — > IR and 
: 57 R d , let q ± := q\ endK ± and ± := </>| en 9 K ±, and set 

MU^^ + T), {0}U:= + ^-), 

[<z]| e := g+n+ + ^n", [0]| e := 0+ • n+ + 0- • n ; 

if e C 9k n 951, we set {</>}| e := 4> + and \q\\ e := g + n + . Finally, we introduce the 
meshsize h : 51 —¥ R, defined by /i(x) = diamK, if a; G k\8k and /i(x) = {/i}, if 

ier. 

Consider the IPDG semi-linear form B(-, ■) : S p xS p —> R, introduced in [HRS05] 
for the solution of the corresponding steady-state problem, defined by 

fl(iu,u):=V / a(w)-\7vdx+ [ (6{a(t,x,h- l \[w]\)Vv} ■ [w] 

(5) „ er ^ ■ /r 

- {«M} • [v] + a[w] ■ [v]) ds, 

where a(w) := a(t,-,\Vw\)Vw, w £ H 1 ^) + S p , for 6 £ {-1,0,1}, with the 
function a : T — >• R + defined piecewise by er| e := C a p 2 /(h\ e ), for some sufficient 
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large constant C a > 0. The corresponding energy norm |||-||| is defined |||iu||| := 

(EkgT II v HIk + Jr^M 2 ^) 172 ' for w G H 1 {n) + S p . The (spatially semidiscrete) 
interior penalty discontinuous Galerkin method (IPDG) for the initial/boundary 
value model problem reads: 

(6) find U : (0, T] ->• S p such that (U t , V) + B(U, V) = (f, V) Vt € (0, T], V € S p . 

3. A POSTERIORI ERROR BOUNDS 

For w G H 1 ^) + S p , and T > 0, we define the norm Hltulll.La^Ti.H'^fi)) : = 
(Jo IIIHH 2 ) ,t>0. We shall derive a posteriori bounds for the error — [7111^2(0 T;ff i(si)). 

3.1. Definition (elliptic reconstruction). Let U be the (semi-discrete) solution 
to the problem ((6]) and fix t G [0,T]. We define the elliptic reconstruction w = 
w(t) G Hq(U) of U to be the solution to the elliptic problem 

(7) (a(w),Vv) =(g,v) Vv € H^Q), 

where g = g(t) is given by g := —AU + f - IT/, with IT : L 2 (Q) S" p is the 
orthogonal L 2 -projection operator onto S' 1 ' and A = ^4(i) : S p —> S p is the discrete 
operator defined by 

(8) for Z e S p , (-AZ, V) = B(Z, V) W G S p . 

The construction of w and that of AZ are both well defined in view of the elliptic 
problem's unique solvability and the Riesz representation, respectively. 

3.2. Remark. The key property of the construction in Definition 13.11 is that U 
is the IPDG solution of an elliptic problem with analytical solution w. Namely, for 
each fixed t g [0, T] it satisfies 

(9) find U E S p such that B(U, V) = (g, V) W G S p . 
We can now decompose the error as follows: 

(10) U — u = p — e, with p := w — u, and e := w — U, 
where w = w(t) denotes the elliptic reconstruction of U = U(t), t G [0,T]. 

3.3. Lemma (differential error relation). Let u, w, U, e, p, e as above. Then, 
for all v G Hq(Q), we have 

(11) (e t ,v) + (a(w) - a(u), Vv) = 0. 

Proof We have 

(e t ,v) + (a{w) - a(u), Vv> = (U t , v) + (a{w), Vv) - (/, v) 

(12) = (U t , v) + (g, v) - (/, v) = (U t ,v) + (-AU, v) - (IT/, v) 
= (U t , n«) + (-AU, Uv) - (/, Uv) = 0, 

using |U), ^} and the properties of the L 2 -projection, respectively. □ 
We consider further the decomposition of U into conforming and non- conforming 
(discontinuous) parts U = U c + U d , where U c G S p nH^(Q) and U d := U~U C G S p . 
Note that there are many ways of performing this decomposition (e.g., by projecting 
U onto the conforming space) whereof the specific nature remains at our disposal 
until further. 

We also use the shorthand notation e c := U c — u and e c :— w — U c ; note that 
e c = p - e c , e = e c + U d and that e c G H^(Q). 
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3.4. Theorem (abstract a posteriori energy-error estimate). With u, U , U d , e, 
and e as defined above, the following error estimate is satishcd: 

, WellUw,™) <Ci|||e|||L2(o,T,^)) + fl~*(||«o - U(0)\\ + \\U d (0)\\) 

+ Cl|||^ HU^O.T.H 1 ^)) + C2||C/t <1 ||L 2 (0 : T : ff- 1 (O)) I 

with Ci := 1 + VZaoT 1 and C 2 := V^a" 1 . 
Proof Set v — e c in (fTT|) . to deduce 

(14) (e$, e c > + (a(C/ c ) - a(u), Ve c ) = -(C/ t d , e c ) + (a(C/ c ) - «(«;), Ve c ). 
Conditions ([3]) and © imply, respectively 

(a(U c ) - a(ji),Ve c ) > a\\ Ve c || 2 , and (a(£/ c ) - a(iu), Ve c ) < a|| Ve c || || Ve c || , 

and the duality pairing (H~ l ,Hl) gives \(Uf,e c )\ < \\U?\\-i\\ Ve c || . Using the last 
3 relations on (|14[) . we deduce 

(15) (e^e c ) + a||Ve c || 2 < (|| U}\\ _ x + a|| Ve c || ) ||Ve c ||, 
which, in turn, implies 

(16) (e^,e c ) + §||Ve c || 2 < ^(ll^ll-i +S||Ve c ||) 2 . 



(17) 



Integrating (|16|) with respect to f between and T, yields 

IKWII 2 + a f ||Ve c || 2 < || e c (0)|| 2 + - f {\\Ut\U + a\\Ve c \\)' 
Jo "Jo 

or 

||Ve c || 2 )' <a"^|| e c (0)|| +^{[ (WtW-i +«||V£ C ||)' 
<a-|||e c (0)||+C 2 ||f/ t d |U 2(O) r,H- 1 (n)) 

+ (Cl - l)lll eC |||L 2 (0,T,ifi(^)) 

noting that |||e c ||| = ||Ve c ||. Using the bounds |||e c ||| < |||e||| + \\\U d \\\, ||e c (0)|| < 
||e(0)|| + ||J7 <i (0)|| on (fl~7|) and the resulting bound on the triangle inequality 

|l|e|||L2(0,T,// 1 (^)) - lll £C HU 2 (o,T,Hi(o)) + \\\U d \\\ L 2(0,T,Hi(n)), 
yields the result. □ 
For the above result to yield a formally a posteriori bound, we need to estimate 
lll e lll-L 2 (o,T,ii' 1 (n)) further. In particular, in view of Remark 13.21 we require an a 
posteriori error bound for the IPDG method for the corresponding elliptic quasi- 
linear problem (j9j. Such a result is available in [HSW] . an instance of which and 
is presented next. 

3.5. Theorem ( HSW]). Let w € Hq(^) be the elliptic reconstruction defined 
in (0) and let W € S p be the solution of (0). Then, for C a > 1 sufficiently large 
the bound 

(18) HI w - IU||| 2 < 8(W, g, SP) := C cst £ (t? 2 + 0(g, W)) , 

holds, with 

rjl = %% + V • a{W)) \\l + -||nr[«(TU)]|| 2 K \ao + C 2 a f\\[W}\\l K , 
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and 

'hi,.,. - w , /' 
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0(g, w DG ) = ^ ( -HI (0 - H) ( 9 + a(W)) \\l + y || (B - U r )[a(W)} ]]dAdQ 



where D denotes a generic identity operator, U denotes the L 2 -projection operator 
onto S p ~ 1 , fir is defined piecewise by firf| e := Tff _1 i>> for al/ elemental faces e cT, 
v e L 2 (fl), where tx p ~ x : L 2 (Q) — > V v -\(e) denotes the L 2 -projection operator of 
the trace on the face e of a function in S^ -1 (with 'Pp^i(e), for e C h the space 
of mapped univariate polynomials of degree at most p — 1 on e), and C cs t > is 
independent of Co-, 0, h and p. 

Also, it is possible to further estimate the terms involving U d , to avoid computing 
U d explicitly. This is done (with, crucially, explicit dependence on p) using the 
following result based on |BE07[ Lemma 3.2]. 

3.6. Lemma. Suppose T does not contain any hanging nodes. Then, for any v € 
S p and any multi-index 7, with \ j\ = 0, 1, there exists a function v c € S p fl Hq(H) 
such that 



(19) ^ \\D^(v-v c )\\i< C 3 \\(-y~ ,u [v]\\ A r, 
with C3 > depending on the maximal angle of T only. 

Proof [BE07, Lemma 3.2] implies that for every k£T there exists an Oswald-type 
operator I 0s : S p ->■ S p n H%(£1), such that 

(20) \\v-I 0s v\\l<C ^IIHIIe. 

for all w € S p , with J^k) := {e e L : e n k ^ 0}. Summing over all the elements 
k 6 T, and observing that the maximal angle and the lack of hanging nodes gives 
an upper bound on the cardinality of J-(k) for all k £ T, we deduce that 

(21) £ii«-*>s<<c£^iimi& 

kGT ecr p 

which shows (fll?)) for I7I = 0. To show (TTI)]) for I7I = 1, we observe that (v — Io s v) G 
S p ; thus, the standard inverse estimate yields: 

(22) £ ||V(« - /os«)|| 2 , < CJ2 £\\v ~ Io s vf K < CJ2 £-\\[v]\\l 

kGT kGT k ecr K 

using the shape regularity of T. Setting v c = IosV, the result follows. □ 

3.7. Remark. The assumptions of Lemma I3TB1 pose the following restrictions on 
the hnite element space S p : the use of quadrilateral elements (as the tensor-product 
nature of the local elemental bases is of crucial importance here), the exclusion of 
hanging nodes and the uniformity of the polynomial degree. If explicit knowledge 
of the polynomial degree p in the a posteriori bounds presented in this work is not 
required, then these restrictions are not needed in view of [KP03[ Lemma 4.1], i.e., 
triangular elements containing hanging nodes can be employed. 

Combining the results of Theorems l3.4l and l3.5[ together with the approximation 
properties described in Lemma I3.6[ we obtain an a posteriori error bound in the 
energy norm for the semi-discrete problem (|6]). 
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3.8. Theorem (energy-norm a posteriori bound) . With the notation of Theorem 
13.41 and the assumptions of Lemma \3.6\ the following error bound holds: 

IINIU 2( o,T,ffi(n)) <Ci f £ 2 (U,g,SP)+a-^\\u ~U(0)\\ 
Jo 

( 23 ) + <r*C 3 || (^) 1 [U(0)]\\r + C4VZ[U]\\l>(o,t;V><t» 

h I 

+ ^ 5 H ("2) 2 I 17 *] IU 2 (0,T;i 2 (r)), 

with C4 := C\\/Cz/C a , C5 := C2CPF and CW > (the Poincare-Friedrichs 
constant), such that < Cpf||u||, for all v G L 2 (il). 

Proof Combining the results from Theorems 13.41 and 13.51 together with the ap- 
proximation properties described in Lemma 13.61 the result follows. □ 
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